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Abstract 

Temporal networks are such networks where nodes and interactions may appear and disappear at various time scales. With 
the evidence of ubiquity of temporal networks in our economy, nature and society, it's urgent and significant to focus on its 
structural controllability as well as the corresponding characteristics, which nowadays is still an untouched topic. We 
develop graphic tools to study the structural controllability as well as its characteristics, identifying the intrinsic mechanism 
of the ability of individuals in controlling a dynamic and large-scale temporal network. Classifying temporal trees of a 
temporal network into different types, we give (both upper and lower) analytical bounds of the controlling centrality, which 
are verified by numerical simulations of both artificial and empirical temporal networks. We find that the positive 
relationship between aggregated degree and controlling centrality as well as the scale-free distribution of node's 
controlling centrality are virtually independent of the time scale and types of datasets, meaning the inherent robustness 
and heterogeneity of the controlling centrality of nodes within temporal networks. 
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Introduction 

The recent outbreak of the A(H7N9) bird flu has caused much 
panic in China, and most of us still remember the financial crisis 
stretching from the USA to the world just a few years ago. These 
two impressive events are typical examples of complex networks in 
our economy, nature and society. Fortunately, considerable efforts 
have been dedicated to discovering the universal principles how 
structural properties of a complex network influence its function- 
alities [1-4]. Not limited to understanding these statistical 
mechanics, another urgent aspect is to improve the capability to 
control such complex networks [5-10], and recent years have 
witnessed the blossoming studies on structural controllability of 
complex networks [11-21]. Classically, a linear time-invariant 
(LTI) dynamical system is controUable if, with a suitable choice of 
inputs, it can be driven from any initial state to any desired final 
state within the finite time [22-24]. Structural controllability of a 
linear time-invariant system, initiated by Lin [25] and further 
developed by other researchers [26-30], assumes free (non-zero) 
parameters of matrices A and B in 

xit) = A'x(t) + Bu(t) (1) 

cannot be known exactly, and may attain some arbitrary but 
frxed values. A directed network, denoted as G{A,E), associated 
with the above LTI system {A ,B) is said to be structurally 
controUable, if ,B) is controUable with the existence of matrices^ 
and B structurally equivalent to A and B, respectively. Noting that 



matrices A and B can be arbitrarily close to A and B when {A ,B) 
is structuraUy controllable, and structural controUabUity is a general 
property in the sense that almost all weight combinations of a 
given network are controllable, except for some pathological cases 
with zero measure that occur when the parameters satisfy certain 
accidental constrains [1 2,25,26] . In the existing literatures [11,12], 
extensive efforts have been focused on the minimum number of 
input signals of such a network. Based on Lin's structural 
controUability theorem [25], Liu et al. [12] stated that the 
minimizing problem can be efficiently solved by finding a 
maximum matching of a directed network, regarding a topolog- 
ically static network as a linear time-invariant system. That is to 
say, a maximum subset of edges such that each node has at most 
one inbound and at most one outbound edge from the matching, 
and the number of nodes without inbound edges from the 
matching is the number of input signals required for maintaining 
structural controUabUity. With the minimum input theorem, many 
contributions to structural controUabUity of complex networks 
have been presented [13-21]. Wang et al. [13] proposed to 
optimize the structural controUabUity by adding finks such that a 
network can be fuUy controlled by a single driving signal. Liu et al. 
[14] further introduced the control centrality to quantify the 
controUability of a single node. Nepusz et al. [15] evaluated the 
controUability properties on the edges of a network. Besides, 
controUing energy [16], effect of correlations on controUabUity 
[18], evolution of controllability [19], controUabUity transition [20] 
and controlling capacity [21], have flourished very recently. 

In our daily life, many networks fundamentaUy involve with 
time. The examples include the information flow through a 
distributed network and the spread of a disease in a population. 
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Table 1. The temporal network in Fig. 1 with the node pairs and active contacts. 





Node Pair(Contact) 


Active Time Points 


Node Pair(Contact) 


Active Time Points 


(A, B) 


[1,2,3,4] 


(B,C) 


[4,6] 


(C, D) 


[2,3] 


(D, E) 


[3,4,5,6] 


(E, F) 


[1,3] 


{B, F) 


[5,6] 


(C, F) 


[4,5,6] 







doi:1 0.1 371 /journal.pone.0094998.t001 



Development of digital technologies and prevalence of electronic 
communication services provide a huge amount of data in large- 
scale networking social systems, including face-to-face conversa- 
tions [31,32], e-mail exchanges and phone calls [33,34] and other 
types of interactions in various online behaviors [35,36]. Such data 
are collectively described as temporal networks at specific time 
scales, where time-stamped events, rather than static ones, are 
edges between pairs of nodes (i.e. individuals) [37]. More and 
more evidences indicate that the temporal features of a network 
significandy affect its topological properties and collective dynamic 
behaviors, such as distance and node centrality [38,39], disease 
contagion and information diffusions [40,41], characterizing 
temporal behaviors and components [42-44] and scrutinizing 
the effects and characteristics within different time resolutions [45- 
47], which are interdependent on the edge activations of temporal 



networks. However, to our best knowledge, a systematic study on 
structural controllability as well as its characteristics of temporal 
networks is still absent. In this paper, similar to the description of a 
static network by a LTI system [12,25], a temporal network is 
associated with a linear time-variant (LTV) system as: 



Mt) = A (t)x(t) + B{t)u{t) 



(2) 



T>NX N 



where A (t)e\ 
matrix of a temporal 
x(t) = ixi(t),X2(t), 



denotes the transpose of the adjacency 
network, i.e.. A' (t) = (A(t))^ , 
captures the time-dependent 

T,NXM 



XN(t)fe\ 



vector of the state variables of nodes, B{t)eM. is the so-called 

input matrix which identifies how external signals are fed into the 




Figure 1. The sequence of graphs representation of the contacts in Table I. In each discrete time point, the network has a different 

formation shown as G\ - ■■ ,G^. 

dol:1 0.1 371/journal.pone.0094998.g001 
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nodes of the network, and u(t) = (ui(t),U2(t), ■ ■ ■ ,UM(t)) e[R is 
the time-dependent input vector imposed by the outside control- 
lers. Meanwhile, by finding and classifying Temporal Trees of a 
temporal network into different types with a combinational 
method of graph theory and matrix algebra, we introduce an 
index as the so-called controlling centrality to quantify the ability 
of a single node in controlling the whole temporal network. With 
analytical and experimental bounds, we point out the indepen- 
dence of the relationship between aggregated degree and 
controlling centrality, as well as the distribution of this centrality, 
over different time scales. Besides, our method reserves as much 
temporal information as possible on structural controllability of 
temporal networks, which may shade new light on the study of 
structural controllability as well as its characteristics without 
wiping out information of the temporal dimension. 

Results 

A temporal network may include a sequence of graphs defined 
at discrete time points. Given a set of N nodes, we denote the 
sequence of graphs as Q = {G' ,t=l,l, - ■ ■ ,T}, where T is the 
sequence length, and G' is a static graph sampled at time point t. 
The adjacency matrix of a temporal network, Q, can be denoted 
by a X time-dependent adjacency matrix A{t\ [=1,2, ■ ■ ■ ,T, 
where a//(0 are the elements of the adjacency matrix of the f''' 
graph, G'. 

For example, a temporal network, Q, with the set of contacts in 
Table 1 can be sampled as a sequence of graphs at time points 
t= 1,2, • • ■ ,6, denoted as 5 = {G^G^G^G^G^G'^} and shown in 
Fig. 1. We illustrate the propagation process taking place on the 



temporal network as shown in Fig. 2. Actually, a message can only 
arrive at nodes B, C and F (dotted nodes in Fig. 2) if its source is 
located on node A, though each node can receive the same 
message if the source is located on node D. This asymmetry (node 
D reaches node A, while not vice versa) mainly due to the direction 
of time evolution, highlights a fundamental gap between static and 
temporal networks. 

2.1 Structurally Controlling Centrality of Temporal 
Networks 

Generally, non-zero entries of a matrix A are free, and A is 
structured if the free entries are (algebraically) independent. Two 
matrices A and A are same structured if their zero entries coincide. 
Matrices A,B,C, - - are independent if all free entries of these 
matrices are (algebraically) independent. In particular, any 
independent matrix must be structured, and any two entries of 
two matrices must be distinct [25,30]. A temporal network is said 
to be structurally controllable at time point to if its associated LTV 
system described by Eq.(2), with a suitable choice of inputs u{t), 
can be driven from any initial state to any desired final state within 
the finite time interval [?o,?i], where the initial and finial states are 
designated at time point ?o and tf {ta<tf <ti), respectively. 

For simplicity, we focus on the case of a single controller and 
reduce the input matrix B(t) in Eq. (2) to the input vector b^"^ with 
only a single non-zero element, and rewrite Eq. (2) as 



x{t) = A!{t)x{t) + b^°^u{t) 



(3) 




Figure 2. The illustration of information propagation on a temporal network, (a), (b), (c) and (d) denote different networks at different time 
points, respectively. Red (gray) time points on edges denote the elapsed time, and the black (dark) time points denote the forthcoming time. 
dol:1 0.1 371/journal.pone.0094998.g002 
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(a) 



(b)\ 




(c) 



(d) 




Figure 3. The illustration of transformation of a temporal network to a static one. (a) Temporal Network with a single controller located on 
node A, (b) The Time-Ordered Graph (TOG), (c) The temporal trees of (a) at time points 1, 2, 3 and 4, (d) the BFS spanning trees of TOG. The red 
(dashed), black (dark) and blue (light) lines stand for the flows of time order, the connection with the single controller and the interactions of 
individuals, respectively. The numbers with parenthesis in (c) denote time stamps. Weights of interactions (the blue ones) are labeled by characters 
«ii,«i2, ■ ■ ■ ,d4i,dn2 in (b), (c) and (d), and without loss of generality, we denote the weight of other edges (the red and black ones) as "1". 
doi:10.1371/journal.pone.0094998.g003 



With non-periodic sampling of Eq. (3), we get its discrete version 
with the recursive relationship for any two neighboring state 
spaces of a temporal network 

x(k + 1 ) = Gi + 1 x(k) + ! u(k),(k = 0, 1 , ■ • • , r - 1 ) (4) 



Define Suio) the structurally controlling centrality of node o in 
a temporal network: 

Sm(o) = rank{ Wc) = rank([GT ■ ■ ■ GjHi , • • ■ ,GtHt- \ ,Ht\) (5) 



where Gk+\=I +Tk+\Ai^^^, Hk+\ = Tk+\U''\ Ai^^^ is the 
transpose of the adjacency matrix of the (k-\- l)th graph, / and 



respectively. Sm(o) is a measure of node o's ability to structurally 
control the network, i.e. the maximum dimension of controllable 
subspace (see Methods), and in this paper, Gk+\ and iZ/c+i are 
structured matrices of size N x N and x 1 , respectively. 

2.2 Graph Characteristics 

Given a temporal network Q{Vg,Eg), where Vg and Eg are the 
collection of nodes and interactions, respectively, we associate G 
with another acyclic digraph N{Q,T). The vertex set of N{Q,T) 
consists of T+ 1 copies, i.e., iiJi, ■ ■ ■ , and ir+i, of each vertex 
ieVg, and T+l copies, i.e., /q,/", • • • , and /f., of the single 
controller denoted as the red ones in Fig. 3 (b). The edge set of 
N(Q,T) consists of three types of edges: (i) the edges connecting 
node i at neighboring time points, ie., it^it + \, t= 1,2, • ■ ■ ,T, for 
each node i^Vg, (ii) the edges it-^Ji + i, where 
i—*jeE(ji, t=\,2, ■ ■ ■ ,T and (iii) the edges connecting the controller 
ie., I°^Oi+ i,/ = 0,l, ■ ■ ■ ,T, where oeVg denotes the directly 
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(a) 



(b) 



(c) 



Figure 4. Three examples of the homogeneously structured temporal trees, (a) Independent trees, (b) and (c) Interdependent trees. For the 
two homogeneously structured trees in (b), there are three same interactions, i.e (B,C,5), (B,D,5) and (B,E,5), but there are only two such interactions, 
i.e (B,C,5) and (B,D,5), for the trees in (c). The trees in (b) and (c) are both interdependent according to our definition. The numbers in parenthesis 
denote active time points of interactions and characters a\,a2,h,h^,h2,c,c\,C2,d,d\ and dj denote the weights of interactions. 
doi:1 0.1 371 /journal.pone.0094998.g004 



controlled node. These aforementioned three types of edges are 
denoted as the red dotted ones, the blue ones and the black ones in 
Fig. 3 (b), respectively. Such interpretation of a temporal network 
is called the Time-Ordered Graph (TOG) model in [39], which 
transforms a temporal network into a larger but more easily 
analyzable static version. For example, we translate the temporal 
network of Fig. 3 (a) to the corresponding time-ordered graph as 
shown in Fig. 3 (b). With the TOG model, we first give the 
definition of input reachability in a temporal network. 

Definition 1. Consider subset = {I^J^, ■ ■ ■ and node 
!V+ie5'2 = {lr+i,27-_|_i, • ■ • ,| oi N(Q,T), which corre- 

spond to node 7° and node ieVg of 5, respectively. If in N{Q,T) 
there exists a path to ij + \ , whose tail I^eSi, then node / of 5 is 
reachable from node 1° at time t, and the set of such reachable 
nodes in Vg is the reachable subset of the input signal /" of Q. 

Proposition 1. The reachability of the input signal of Q is 
equivalent to the reachability of subset Si, i.e. the (t+l)''' row of 
{T+l — t)'^' power of adjacency matrix of N{Q,T), and the 
controlled rows of dynamic communicability matrices of Q starting 
at different time points t, denoted as {2'}oV) where 

r=l,2,--,r. 

Proof. Denote partitioned matrix An^q^t) (size (iKgl + l) 
x(|Kg| + l)) as the adjacency matrix oi N{Q,T), and for each 
block B(jj^ (size (T+ 1) x (T+ 1)) of matrix Afj^g^j), if there's a 
directed edge ii—^jt+i in N{Q,T), where t= 1,2, ■ • • ,T, then we 
have {B^ij)},,^^ t^O and {AN(g,T)} nT+\)+ij(T+\)+i+\^'^- RecaU 
the dynamic communicability matrix [40] to quantify how 
effectively a node can broadcast and receive messages in a 
temporal network, defined as: 



Here, matrix Ai is the adjacency matrix of the graph, and 
0<fl<l/p (p denotes the maximum spectral radius of matrices). 
Similarly, we define the communicability matrix starting at 
different time points to quantify the reachability of the controller, 
written as: 



Q, = {I'' +a,A*,)(r +a, + lA*+i )■■■ {r+arAlj,) 



(7) 



where A* - 



is the adjacency matrix of the / 



Note that a non-zero element {ij) of 

k 



0 (M"))' 
.Oatxi At 

graph with a single controller 1° located on node o, and 
0 Oixw 

a product of matrices, such as {A^ , is the reachability from node / 
to node j if {(^)^}y7^0, and the length of paths in graph N(Q,T) 
is never more than 7" + 1. Therefore, the reachability of node 1° in 
Si = {/°,/°, ■ ■ ■ ,7°} is the {t+ l)'* row of {T+ \-t)"' power of 
AN(g,T), i.e., {{AN(g.T)f^^''}i+iy, where t = 0,\,---,T. For 
each column of matrix Wc, we have Gt ■ ■ ■ GjHi = 
[(Gt ■ ■ ■ G2)]Hi = [G2 • • • Cfj-]Hi = [iI+A2) ■ ■ ■ (I+At)]Hu 
and with the definition of matrix Qt, we know that {Qi}i j 
describes the reachability from node / to node j. Therefore, the 
rechabUity of controller 1° at time t is equivalent to the controlled 
row, i.e. the o"" row, denoted as {Qt} of matrix Qt. 

With Proposition 1, we rewrite matrix Wc in the form of 
reachability as: 



--(I + aAi)(I + aA2)---{I + aAT) 



(6) 



w*=[({QiU) x{e2U),---,({erU): 



(8) 
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Figure 5. Controlling centrality of artificial networlts. (a), (b), (c) and (d) denote network with 40, 60, 80 and 100 nodes, respectively. For each 
of the four networks, we randomly generate an interaction between a pair of nodes with probability 0.002, and repeat it for all the N{N — l)/2 pairs 
of nodes at a specified time point, repeat this process for 100 rounds at 100 different time points, i.e. t= 1,2, ■ ■ ■ ,100. The value of controlling 
centrality, denoted as 'Calculated', is straightly calculated by the computation of matrix IV^- in Eq. (19), and the upper and lower bounds, denoted as 
'Upper Bound' and 'Lower Bound', respectively, are given by the analytical results in Eq. (16). 
doi:1 0.1 371 /journal. pone.0094998.g005 



where v denotes the reachability of the controller at time 

point t, and we have rank(W*) = rank(Wc). As shown in Fig. 3, 
we easily get {gi}i_v = [0, 1 +a2iC3i^42, a2i, a2iC3i+^4i, fl2i'^4il, 
{e2}i,v = [0,lA64i',0], {gsli.v =[0,l,OMuO], and {Q4}i, = 
[0,1,0,0,0]. According to Proposition 1, = [{{Qi}iy)' , 

({e2}i,v)',({e3}i,v)', 



({e4},,v)'] = 



/ 


0 


0 


0 






1 +a2iCT,\bn2 


1 


1 


1 




aix 


0 


0 


0 




fl2lC31 +^41 


641 


^41 


0 


\ 




0 


0 


0/ 



Definition 2. A temporal tree, denoted as TT,, of a temporal 
network Q{ Vg,Eg) is a Breadth-First Search (BFS) spanning tree, 
denoted as STi, of its corresponding static network N{Q,T) (TOG 
model) rooted at node /feSi = ■ ■ ■ ,/?■}. 

Remark. The Breadth-First Search (BFS) is a classical 
strategy for searching nodes in graph theory, and a BFS spanning 
tree contains all the nodes and edges when the BFS strategy is 
applied at some node. A distinctive property of N{Q,T) is that 
there's no cycles in it, and each path's length is no more than 
7^ + 1, so it's easy to apply the BFS strategy to find trees rooted at 



some designated nodes in N{Q,T). Obviously, the one-one 
mapping between a temporal tree of a temporal network and a 
BFS spanning tree of the TOG is guaranteed by the one-one 
mapping between Q{Vg,Eg) and N(Q,T). For the temporal 
network in Fig. 3 (a), each of the three temporal trees, as shown 
in Fig. 3 (c), of this temporal network exists a unique 
corresponding BFS spanning tree, as shown in Fig. 3 (d). 

Proposition 2. Denote RtTi^RtTi^' ■ ■ ^ and Rjtt^ 
|[j(K(; + l)xi jj^^ reachability vector of each temporal tree from 
the controller 7°, and matrix = [RtTiRtT2' ' ' ^tTt\ 

glRt^-e + Dxr^ we have rank{W'^) = rank{W*). 

Proof. With Proposition 1 and Definition 2, we know there's a 
temporal tree 7"7", of each ST^ in TOG, and each ST^ is a leading 
tree when compared with ST,j^\ (refer to the definition of BFS 
spanning tree with the TOG model). Therefore, each temporal 
tree TTt is a leading tree when compared with TT,j^\ . Two 
strategies are adopted to yield a leading temporal tree: i) Adding 
new nodes into TT,, i.e., we have | f^rr, I > I ^^ttz+i I, ii) Adding 
new paths to the existing nodes, i.e., we have | Vjt, I = I Vtti+ 1 1 ■ In 
the case of strategy i), if there's only one temporal tree, we 
obviously have rank(W^) = rank(RTT) = rank(W*)=l; if the 
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Figure 6. The gap of upper and lower bounds of controlling centrality. (a) HT09 (b) SG-lnfectious (c) Fudan WIFI. For the dataset of 'HT09', 
two temporal networks are generated: I) a temporal network (113 nodes and 9865 Interactions) with all nodes and Interactions within record of 
dataset, denoted as 'all range', II) a temporal network (73 nodes and 3679 Interactions) with nodes and Interactions after removing the most powerful 
nodes (nodes with the largest controlling centrality) In the temporal network of I), denoted as 'removed'. For the dataset of 'SG-lnfectlous', three 
temporal networks are generated: I) a temporal network (1321 nodes and 20343 Interactions) with nodes and Interactions recorded In the first week, 
denoted as 'Week 1', II) a temporal network (868 nodes and 13401 Interactions) with nodes and Interactions recorded In the second week, denoted as 
'Week 2', III) a temporal network (2189 nodes and 33744 Interactions) with nodes and Interactions recorded In the first two weeks, denoted as 'Week 
1 &2'. For the dataset of 'Fudan WIFI', three temporal networks are generated: I) a temporal network (1 1 20 nodes and 1 2833 Interactions) with nodes 
and Interactions recorded In the first day, denoted as 'Day 1', II) a temporal network (2250 nodes and 25772 Interactions) with nodes and Interactions 
recorded In the second day, denoted as 'Day 2', III) a temporal network (1838 nodes and 27810 Interactions) with nodes and Interactions recorded at 
Access Point No. 71 3, denoted as '71 3 point'. The upper and lower bounds of the controlling centrality are given by analytical results In the main text, 
and the gap Is given by the absolute value of the difference of the upper and lower bounds. The aggregated degree of a node Is the number of 
neighbored nodes whom It Interacts within the corresponding temporal network. All the gaps are minor when compared with the sizes of these 
temporal networks. 

dol:10.1371/journal.pone.0094998.g006 



number of temporal trees is n, and rank(W*) = rank(W^) = n, 
then when the number of temporal trees is « + 1 , we have 

/Oix„ 0 \ 

rankiW'^) = ranki W*^„ 0)„xi \=rank{W*) = 

\%Vg\-n)xn (t)(|Ks|-«)xl / 

« + 1 , where ({) denotes a nonzero vector. In the case of strategy 
ii), each new interaction in leading tree TT,, which isn't included 
in temporal tree TT, + 1 , contributes to new paths to the existing 
nodes. By some linear superposition of columns of matrix W* and 
W^, we find there's no impact on the maximum rank of matrix 
W* if we cut down and drop those "old" interactions, which 
means we only need to take the leading temporal 
tree, i.e. TT^, into consideration. Therefore, we have 
rank{ W'^) = rank{Px W'^Tx)^ rank( W*) = rankiPj W*T2), 
where Pi, Pj, T\ and are properly defined linear transforma- 
tion matrices. 

For example, according to Definition 2, the reachability of 
temporal tree TTi of Fig. 3 (c) is RtTi = [0,l,fl2i,a2iC3i,a2i'inl ■ 
Similarly, we have -Rt-tt, = [0,1,0, C4i,0] , /^ttj = [0,1,0, C4i,0] and 
iJj,j,^ = [0,1,0,0,0]' for temporal trees TTj, TTj, and TT4, 
respectively. Therefore, we easUy reach P] W'^T] =Pi [RtTi ,RtT2> 



Rtt,,Rtt.]T\ 



I 0 
0 

021 
fl2lC31 

\a2in'4i 

/ 

1 + 
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0 
0 

641 
0 



Pi 
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0 

021 f 31*42 
fl21 

ai\c^\ +^41 

«2lfi?41 



0 
0 
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0 
0 
0 

b^\ 
0 



/ 0 
1 

021 

fl2lC31 

0\ 
1 
0 
0 

0/ 
0 

1 

0 

&41 

0 

o\ 
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0 

0/ 



0 
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0 

/341 

0 



0 
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0 

641 

0 



1 

0 
0 

0/ 



and P2WT2-- 



0 
1 
0 

641 
0 



1 

0 
0 

0/ 



T2 = 



.Obviously, rank(W'^) = rank{Pi W'^Ti) 



--rank{W*) = rank(P\ W*Ti), which is consistent with Proposition 



2. 



PLOS ONE I www.plosone.org 



7 



April 2014 I Volume 9 | Issue 4 | e94998 



SC and CC of Temporal Networks 





120 n 




100- 




80- 




60- 




40- 




20- 




0- 







^ Lower(all range) 
Upper(all range) 

□ Lower{removed) 
Upper(removed) 



^^^^ 



0) 


100-, 


O 


80- 


D) 
C 


60- 




40- 




20- 


n 




o 


0- 


O 





— I— 

20 



^ Lower(week 1 ) 
Upper(week 1 

□ Lower(week 2) 
Upper(week 2) 

O Lower(week 1&2) 

^ Upper(week 1&2) 



D) 500 
% 400. 
> 300 ■ 
200 
100 
0 



^ Lower(day 1 ) 
Upper{day 1) 

□ Lower(day 2) 
Upper{day 2} 

O Lower(713 point) 

^ Upper(713 point) 



— I— 

10 



—1 — 
40 



— I— 

60 



—1 — 

80 




— I— 

20 



— I— 

30 



40 



100 200 300 

Aggregated Degree 



400 



100 



500 



(a) 



(b) 



50 



(C) 



Figure 7. The statistical relationship between node's aggregated degree and the average controlling centrality. (a) HT09 (b) SG- 

Infectious (c) FudanWIFI. All the temporal networks are the same as those in Fig. 6. Each point in this figure is an average controlling centrality of 
nodes with the same aggregated degree, and there's a positive relationship between the aggregated degree and its controlling centrality, even with 
some structural destructions or time evolutions. 
doi:1 0.1 371 /journal.pone.0094998.g007 



Definition 3. Temporal trees TT\,TT2,--- are homoge- 
neously structured if their corresponding adjacency matrices, 
denoted as Att\, ^tTi.' • are same structured. Otherwise, they 
are heterogeneously structured. 

We rewrite matrix as: 



(9) 



In Eq. (9), matrix W'' = ( Oj-z. x i( )' ) of size (| Vg\ + 1) x 
denotes the part of heterogeneously structured trees, and matrix 
W'^ = ( Ojs X ) of size (I Vg\ + 1) X T'^ denotes the part 

of homogeneously structured trees, respectively. Obviously, 

2.2.1 Heterogeneously Structured Trees 

Definition 4. If heterogeneous trees TTd],TTo2 , ■ ■ ■ consist 
of same nodes, i.e., Vttdi = ^tTdi = " then they are called 
heterogeneous trees with same nodes. Otherwise they are 
heterogeneous trees with different nodes. 

To determine the rank of matrix , we rewrite it as: 



^2 



(10) 



In Eq. (10), each W^^ 
{\Vg\ + l)xTg^ denotes a collection 
trees with same nodes ( V^yo. ^ VfyD, 



(OrlxiCW'^)') ,/=l,2,---,;c, of size 
of heterogeneous 
for V/t^/'), and 



f''t, = (07-f^„xi(»'t,)')' = [i^?r,^^?r.^---.^?V ] of size 
(I Fgl + 1) X r^^j denotes heterogeneous trees with different 

nodes. Y.ilih, = T''. 

Case 1. Heterogeneous trees with same nodes. 

Proposition 3. Given matrix as a collection of 

heterogeneous trees with same nodes, we have 
rank(W§') = min(\V^yD-lT§), /= 1,2, • • • where \V^b>\ de- 

notes the number of nodes in matrix W^^ . 

Proof. According to the definition of heterogeneous trees with 
same nodes, these trees always have the same reachability with 
dilferent paths to reach the same node, which means for each 
heterogeneously structured temporal tree with same nodes, there 
exists at least one independent parameter (interaction). When 
r^ = l, raMA:(frJ^*W)=l. When r^ = w<|F>^'|- I, we get a 
triangular matrix with its diagonal elements non-zeros by some 
linear transformations. Therefore, we have rank( Wg^ ) = Tg^ = n. 

Similarly, when | Vwf \<T5=n, we get rank( W^' ) = \ K^b* | . 

' ^/ 

In short, we reach rank{W^ ) = min(\ F[yD« \,T^). 

Case 2. Heterogeneous trees with different nodes. 

Proposition 4. Given matrix as heterogeneous trees 

with different nodes, we have rank(W^'^^^) = 7)^^^ . 

Proof. When Tf^^, = 1, we easily have rank( W^^''^^ )= 1. If 
TS =m, and rawA:( ) = M, then when TS =m+1 it's 
equivalent to add a tree with different nodes into matrix Wg^ to 
get matrix ■ Therefore, there always exists 

at least one new nonzero entry with its column index 
n+l and row index r(n<r<\Vg\ + l) in matrix 



W. 



St + i 



and rankiwf"^'^) = rank {[W^'^"' R?^] ) 
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Figure 8. The specific relationsKiip between node's aggregated degree and controlling centrality. (a) and (b) Temporal networks 
generated by the dataset of 'SG-lnfectious' (c) and (d) Temporal networks generated by the dataset of 'Fudan WIFI'. Although big nodes (node with 
larger aggregated degree) tend to own larger controlling centralities, there exist many nodes with larger (smaller) aggregated degree but smaller 
(larger) controlling centrality, such as circled points in (a), (b) and (d). 
doi:1 0.1 371 /journal.pone.0094998.g008 



/Oix„ 0 \ 

rank I ^nxn (t)«xi =« + l, where (J) denotes a 

\%Vg\-n)y.n (t){| Kgl -«) x 1 / 

nonzero vector. That means for any Ts^_f_^, we have 
rankiWS ) =TS . Note that if we cannot find such a nonzero 
entry, we claim that this new tree must have a collection of nodes 
coincident to some other tree, which is not allowed in this case. 

Theorem 1. Given matrices Wg^ , /= 1,2, • • • 1, as the 
heterogeneously structured trees and Sm{o) as the maximum- 
structurally controllable subspace of heterogeneously structured 
trees, we have 



ma^+l{rank(W|*)}<S•j:,^^ 

k+i (11) 
= rank( JF^) < ^ rank( wf ) 

1=1 

Proof. Firstly, we prove the left part of inequality (11), i.e. 
rank(W'^)>max{^! {rank(Wsl')}. Compared with the 
trees, denote as TT],TT2,---, TTtP , in matrix WS 
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Figure 9. The distribution of node's controlling centrality. (a) Temporal networks generated by the dataset of 'SG-lnfectious' (b) Temporal 
networks generated by the dataset of 'Fudan WIFI'. For each dataset, three different temporal networks are generated within different time scales, 
denoted as 'Week V, 'Week 2' and 'Week 1&2' for SG-lnfectious and 'Day 1', 'Day 2' and '713 point' for Fudan WIFI, respectively. 
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^^St 1 =[^Tr,' ^77-2' ^TT r, D' ''"'^^^ ™ matrices 

Wg^ ,1 = \,2, ■ ■ ■ ,k, have difTerent nodes, i.e., 

Vwi; # Vtt, ¥=Vtt2¥= ■■■ ¥= Vjt^^ , and Vw^' 7^ V^^' 

for V/t^/'. Therefore, rank(W'^) = max'^:! {rank{W^I)} 
when there exists a matrix consists of all nodes, and it has the 
maximum rank. For the riglit part, i.e. 
rank{W^)< X^f.V rankiWg^ ), we reach the equality when 
matrbc is written as: 

1 1 ... 1 



1 



0„ 



0. 



•^k -I 



W* 



where row vector li x , i.e, the 2* row of matrix W'^ , denotes 
node o, and matrices W*^ ,W*'^ ,■ ■ ■ ,W*^^^ denote the other 
part of these trees. This means there's no intersection of nodes 
between any two matrices of Wg^ ,/ = 1,2, ■••,/: + 1, except node 
o, i.e., I V^B- n V^yD'^ I 



: 1 for Ml^l . In this case, each matrix 



1^1* contributes d = rank(W§^) = T^^ to rank{W'^). Therefore, 
rank(W^)= ^*+/ rank{W^*). 

2.2.2 Homogeneously Structured Trees 

Definition 5. Consider homogeneously structured trees 
TTsi ,TTs2 , ■ ■ ■■ If their corresponding adjacency matrices 
A-TTsi , AjTs^ 5 ■ ■ ■ are independent, then they are called inde- 
pendent trees. Otherwise they are interdependent trees. 



We rewrite matrix as: 



(12) 



and each iy)^,(w= 1,2, ■■ • ,q), denote a collection of homoge- 
neously structured trees (K^/s ^ Vyy^, for Vm ), which is 
written as: 



K- 
/o, 



In 



Eq. 



Sm,p{m) 



Sm,p(m)+l 
'^'^ Sm,p(m)+\ 



(13) 



13), 



each 



W: 



Or? 



w = 1,2, • • • ,p{in), of size (| Kg| 1) x Ti_,n' denotes a collection 
of interdependent trees with same interactions {Im,w¥^Im,w' for 
Vu'T^iv , where /,„,„. denotes the collection of same interactions in 



matrix W^^X and W'X,X'")+l^ 



size (|Fc| + l)x r; 



denotes independent trees. For homo- 
geneously structured trees, we have Xim = i X^hS'i^' '^s„w~'^^ 
and I f^n^ I = I I'^f^'i' I, w'= 1,2, • • • ,/7(ni)+ 1, where l^w^} and 
I Vif rs' I denote the number of nodes in matrices and w 
respectively. 
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Table 2. Notations in the paper. 




Notations 


Description 


G' 


the t'^' formation of temporal network Q(Vg,Eg) 


V and 1 V\ 


the set of nodes and the cardinality of set V 


A, 


the adjacency matrix of graph 


(A) 


the transpose of adjacency matrix A 


{Af 


the A:'^ power of adjacency matrix A 




an element of matrix A with position / (row Index) and / (column index} 




the i'^ row of matrix A 


r 


the controller located on node o of temporal network Q{ Vg,Eg) 


a 


dynamic communicability matrix of temporal network Q{Vg,Eg) at time t 


w 


reachability matrix of input signal within the temporal network Q{Vg,Eg) 


Rtt 


reachability vector of input signal within a temporal tree TT 




reachability vector of input signal within heterogeneously structured 


temporal tree TT 




reachability vector of input signal within homogeneously structured 


temporal tree TT 


W" 


reachability matrix of input signal within temporal trees extracted from 


temporal network Q{Vg,Eg) 


w" 


reachability matrix of input signal within heterogeneously structured 


temporal trees 




reachability matrix of input signal within homogeneously structured 


temporal trees 




the maximum controlled subspace of temporal network Q(Vg,Eg) 


with single controller located on o 




the maximum controlled subspace of heterogeneously structured 


temporal trees with single controller located on o 




the maximum controlled subspace of homogeneously structured 


temporal trees with single controller located on o 


doi:1 0.1 371 /journal.pone.0094998.t002 



Case 1. Independent trees. 

Proposition 5. Given matrix Ws*p(,n)+l as indepen- 
dent trees, we liave rank{Wl'^j^^^^i) = min{\Vu,s\,Tl_p^^^^i), 
where | Viys\ denotes the number of nodes in matrix ^s,„ p(m) + \- 
Proof. According to the definition of independent matrices, 
the matrix having the reachability vectors of independent trees 
from the controller i.e. [J^TTsr^TTsi' ' ' ^ structured 
matrix. For such a structured matrix, we can always find a square 
sub-matrix of size min(\ f^»'il.7'|„„;,(„,)+i) x mm{\ ^wslTl„p(m)+il 
whose elements are all non-zero. Therefore, it's obvious that 
rank( W^^ + 1 ) = V^s \ , + 1 ) . 

An illustrative example is given with Fig. 4 (a). The 
corresponding matrix ^s„ p(m) + l written as: 

0 1 fl] a]hi a\Ci aidi 

0 1 (72 ^^2^2 O2C2 a2'4 



S„,.p{m)+\ - 



whose rank is 



2, i.e., '•««^(W^i:,,(,„, + i) = 7'|„„„„,) + i=2<|F^.|=5. More gen- 



eraUy, if r|, ,^(,„)+r 



: w > I Vws 1 = 5, matrix Wf 



'^'^ S,„,p{ni)+\ ■ 



rank(WH^^„,y 



( 0 


0 




0 


0 \ 


1 


1 




1 


1 


a\ 


fl2 




Cln 


Cln 


a\h\ 


(32^2 


■ • an 


-\bn-\ 


Unbn 


a\C\ 


132 C2 


■ ■ an 


-\Cn-\ 


ClfjCfi 


a\d\ 


02 


Cln 


-\dn-\ 




\ 








n n ^ 




= 5<r| 


,p{m)+\ 


= n. 





Case 2. hiterdependent trees. 
Proposition 6. Given matrix 



and 



as a collection of 



interdependent trees, we have rank{Wf^^^) = min(\Vw^} — 
\hn.w\,T^^^K)X'^v=\,2,---,p{m)), where \ Vfys^\ denotes the num- 
ber of nodes, and |/n7,iv| is the number of same interactions in 

Proof. Without loss of generality, we firstly prove the case of 
two trees as shown in Fig. 4 (b). Here \Im,w\ =3, i.e. interaction 
(B,C,5),{B,D,5) and (B,F,5). The corresponding matrix 

0 1 fli aih flic a]d^ 
0 1 02 uic 
0 1 fli aib aic aid^ 
0 -02/01 0 0 0 0 



, and it's obvious that the 
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Table 3. Characteristics of the three empirical datasets. 






HT09 


SG-lnfectious 


FudanWIFI 


Area 


Conference 


Mesume 


Campus 


Technology 


RFID 


RFID 


WiFi 


Collection 


3 days 


62 days 


84 days 


duration 


Number of 


113 


10970 


17897 


individuals 


Number of 


9865 


198198 


884800 


contacts 


Spatial 


<2 


<2 


<8 


resolution{meters) 


Types of 


Strangers 


Acquaintances 


Acquaintances 


contacts 


with repeat 


without repeat 


with repeat 


doi:l 0.1 371 /journa!.pone.0094998.t003 



dependence of elements in matrix is caused by tlie interdependent 
of trees in some interactions. Thus rank( W§* ^^,) = 
2=\Vws\-\I,. 



■ s,„, 



More generally, when extending to 



the case of n trees, Wg - 



( ^ 


0 


1 


1 


fll 




aib 




a\c 




aid 


ajd 



0 \ 

1 

a J} 

ClfiC 

cif,—\d cif,d 

I 



0 

1 

a„-\b 
an-\c 



( ° 


0 


0 


0 \ 


1 


-02/01 ■■ 


— o„-\/ai 


-a„/ai 


fll 


0 


0 


0 


a\b 


0 


0 


0 


a\c 


0 


0 


0 


aid 


0 


0 


0 


\ 






/ 



\VwS,\- 



\Im,n \ < 2"? „ = 



-n. Similarly, for the trees 



in Fig. 4 (c). 



fll 



oib a\c o\di ^ 
oib a\c o\d2 



1 Oi 

-02/01 0 



oib oic aid\ ^ 

0 0 (£/2-fi?l)fl2 



and 



ranki WlJ = 2 = r|__ „, < | V^s^ \-\I„, 

tending to the 

/ 0 
1 



oib 

OiC 



0 
1 

02 
02 ft 

fl2C 



0\di 02d2 



= 3. Similarly, when ex- 
case of n trees, 

0 0 \ 

1 1 

fl/7 — 1 fl/7 

a„_\b a„b 

fl„_lC fl„C 



v 



/ 0 


0 


0 


0 \ 


1 


-02/01 


-fl„_i/fli 


-fln/fll 


fll 


0 


0 


0 


Oib 


0 


0 


0 


OiC 


0 


0 


0 


oidi 


(d2 — di )a2 


■ (dn-i—di)o„_i 


(d„-di)a„ 


\ 






J 



and W'l;,,, = 3 = \V^s\- 1/„,, J < Tl„ = n. 

Theorem 2. Given matrices Wf w= 1,2, •• ■ ,p(m) + 1, as 
homogeneously structured trees, we have 



p(m) 



rank(W^) = min{min[^^ (rank{Wg^^ i^,)), 
11- =1 

"ifl4=i{| V^s\- |/„Mv|}l +rank(W|l^^,„^^l),\ V^s\} 



(14) 



where \Vw^s \ is the number of nodes, and \Im,K\ is the 
number of same interactions in U'= 1,2, ■■ • ,p(m). 

Proof. The outsider function min{} ensures that the rank of 
matrix W^^ never exceeds the number of independent rows, i.e., 
the number of nodes in matrix fF^. Next we focus on the number 
of independent columns. From the proof of Proposition 5, we know 
there always exists a structured square matrix of size 
mm(\Vfvs\,Ts,„,p(,„)+i)xmin(\Vws\,Ts,„,p(,„)+i) in matrix 



so there always exists mifi(\Viys\,Ts,„,p(m)+l) inde- 
pendent columns compared with interdependent matrix 



S„,p(m)+\ 



W"' 



which means matrix 



S„,p{m)+\ 



always contributes ronkiWf^ to matrix W,^, 



i.e., >'ank(Wf^ in Eq. (14). Now we focus on the part of 

'«'4E:i'i ironk(WlJ),mo4'::\{\ V^s\- 1/„,,,,|}1, which 
deals with the rank of all interdependent trees, i.e. the 
rank of matrix W'd. Without loss of generality, for 
trees shown in Fig. 4 (b) and (c), we have 
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,and rankiWl^ i Wl^ 2) = ^<'''^i^L,i) + ''^(^s„,2) = ^- More 
generally, when extending to the case of n trees, we similarly have 



/ 0 

1 


0 

1 




0 0 

1 1 


0 

1 


0 ■■• 0 

1 ... 1 


1 




"2 






fi 


"2 - ■ 




aib 


a2h 




■ <l.-l* 0.6 


aib 
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arc 




- an-\c a„c 


a\ci 






\axd 


aid 




• a„-id a„d 


a\d\ 


a^di ■■• a„-idn- 


1 a,4, / 






0 


0 


0 


0 


0 


1 --- 




- 




0 










0 


0 


0 


0 


0 






0 


0 


0 


0 


0 






0 


0 


ai(Ci - 


-c) a2(C2-c) ■■■ 








0 


0 


ai{d\ - 


d) aiidi — d) ■■■ 


On-lidn- 



, and rank{ W^^ ^ ^s„,2^ = ^ < rcmk( + rcmk( W^^ j) = 5- 
When E^J?i('-««^(f^llJ)^'«<*=Ul^'»'^|-|^.,«'l). it's easy to 
verify that rank{ Wjl) = Ylw=\ i^'^^i^sl, »))• ^ ""^^ of matrix 
W''' is {rank{Wl J\ma^^:i\{\ - 1/„,,,|}]. 

With Eq. (12) and Theorem 2, we directly give the following 
Lemma 1 for homogeneously structured trees. 

Lemma 1. Given matrices W^, m=\,2, - ■ ■ ,q, as collections 
of homogeneously structured trees and 5^(0) the maximum- 
structurally controllable subspace of homogeneously structured 
trees, we have 

'I 

maxl ^ 1 rank{ W^) < Sf^^^^ = rank{ W^) < ''""^C ^m) (15) 

m=l 

With Theorem 1, Theorem 2 and Lemma 1 above, we straightly get 
Theorem 3: 

Theorem 3. Given SM{q) arid Sm^o) ^ the maximum 
controlled subspace of heterogeneously structured and homoge- 
neously structured temporal trees in Eq. (1 1) and (15), respectively, 

we have 

wax(5^(„),5f,(„)) < 5m(o) < + (16) 



2.3 Numerical Simulations 

We firstly verify the feasibility and reliability of Theorem 3. As 
shown in Fig. 5, four different networks with sizes of 40, 60, 80 and 
100 are studied, respectively. For each of the four networks, we 
randomly generate an interaction between a pair of nodes with 
probability 0.002, and repeat it for all the N{N —\)/2 pairs of 
nodes at a specified time point. Repeat this process for 1 00 rounds 
at 100 different time points, i.e. t= 1,2, • • • ,100. As shown in Fig. 5, 
aU the calculated values of controlling centrality of the four 



networks (denoted as 'Calculated') are between the upper and 
lower bounds (denoted as 'Upper Bound' and 'Lower Bound', 
respectively) given by our analytical results in Eq. (16). Besides, the 
gaps (numerical calculations) between upper and lower bounds are 
very minor in these artificial networks. 

We further investigate three empirical datasets, i.e., 'HT09', 
'SG-Infectious' and 'Fudan WIFE (Details of the datasets see 
Methods) [31,35,36,43]. For the dataset of 'HT09', two temporal 
networks are generated: i) a temporal network (113 nodes and 
9865 interactions) with all nodes and interactions within record of 
dataset, denoted as 'all range', ii) a temporal network (73 nodes 
and 3679 interactions) with nodes and interactions after removing 
the most powerfiil nodes (nodes with the largest controlling 
centrality) in the temporal network of i), denoted as 'removed'. For 
the dataset of 'SG-Infectious', three temporal networks are 
generated: i) a temporal network (1321 nodes and 20343 
interactions) with nodes and interactions recorded in the first 
week, denoted as 'Week 1', ii) a temporal network (868 nodes and 
13401 interactions) with nodes and interactions recorded in the 
second week, denoted as 'Week 2', iii) a temporal network (2189 
nodes and 33744 interactions) with nodes and interactions 
recorded in the first two weeks, denoted as 'Week 1&2'. For the 
dataset of 'Fudan WIFE, three temporal networks are generated: i) 
a temporal network (1120 nodes and 12833 interactions) with 
nodes and interactions recorded in the first day, denoted as 'Day 
r, H) a temporal network (2250 nodes and 25772 interactions) 
with nodes and interactions recorded in the second day, denoted 
as 'Day 2', iii) a temporal network (1838 nodes and 27810 
interactions) with nodes and interactions recorded at Acc:ess Point 
No. 713, denoted as '713 point'. With these three types of eight 
temporal networks, we calculate their upper and lower bounds of 
controlling centrality given by our analytical results (it's difficult to 
directly calculate the rank of matrix Wc for large-scale networks). 
The aggregated degree of a node in Figs. 6 and 7 is the number of 
neighbored nodes whom it interacts within the corresponding 
temporal network. As shown in Fig. 6, although the sizes of these 
networks range from 73 to 2250, the gaps of the upper and lower 
bounds remain very tiny, indicating the feasibility and reliability of 
Eq. (16) in both artificial (refer to Fig. 5) and empirical networks. 
Fig. 7 shows us the positive relationship between the aggregated 
degree and controlling centrality of nodes. When removing the 
most powerful nodes (nodes with the largest controlling centrality), 
as shown in Fig. 7 (a), and considering temporal networks with 
different time scales and types, as shown in Fig. 7 (b) and (c), the 
observed positive relationship remains unchanged. This indicates 
the robustness of this relationship of temporal network, regardless 
of the structural destructions or time evolutions of the network. 
Further more. Fig. 8 reveals some nodes with rather larger 
(smaller) controlling centrality but smaller (larger) aggregated 
degree, which suggests that although there's a positive relationship 
between aggregated degree and controlling centrality, controlling 
centrality is a measurement inherently different from the 
aggregated degree. 

Besides, Fig. 9 focuses on the datasets of 'SG-Infectious' and 
'Fudan WIFI' to visualize the distribution of controlling centrality 
of different temporal networks. The scale-free distribution of 
node's controlling centrality is virtually independent of the time 
period and network scale, which is similar to the distribution of 
node's activity potential [47] . However, these two studied datasets 
are inherently different. The dataset of 'SG-Infectious' collected 
the attendee's temporal activity information during an exhibition, 
where the attendee generally do not appear again after the visit. 
Therefore, the interactions among nodes in the temporal networks 
generated from 'SG-Infectious' present more randomness than 
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those of 'Fudan WIFI', while the latter presents weekly rhythm of 
the scheduled campus activities in a university. 

Discussion 

Many problems on nc'tworks involving time are raised by some 
common themes, especially on communication in distributed 
networks, epidemiology and scheduled transportation networks. In 
some earlier literatures, authors studied a model with each edge of 
a graph associating with a single active time point (or equivalendy 
a single starting and ending time points). So each edge has a 
reaction time, i.e. the delay, to transmit an information to the 
other end of the edge. This simple model had raised a number of 
interesting open questions about the basic properties of the 
original graph. However, such a simplified model is far enough for 
the cases in our information society, where relationships are 
varying, they are strengthened or weakened, even disappeared or 
created, and the exchange of information happens frequently, i.e. 
a pair of node exchanges information for far more than just once. 
Although, with the record of temporal networks being available by 
digital technologies, many attentions have been attracted to this 
area [37], litde work has been carried out on the structural 
controllability. 

In this paper, we propose a framework from graphic 
perspective to address the structural controllability of temporal 
networks, especially focusing on the ability of a single node to 
control the whole network (controlling centrality), which allows us 
analyzing large-scale networks more convenient and efficient. 
Noting that the single node here does not necessarily to be driven 
node as the one in the seminal paper of Liu et al. [12], it is 
randomly chosen from the whole network and we mainly focus 
on its controlling centrality - a measurement of its importance 
from the persjjective of control theory. Although there's a positive 
relationship between controlling centrality and aggregated 
degree, these two centralities are obviously not equivalent in 
neither definition nor methodology. Frankly, more steps can be 
taken on the structural controllability of temporal networks in the 
near future. For example, one of opening problems, untouched in 
this paper and waiting for endeavor studies and explorations 
elsewhere, is the multi-inputs case. Whether or not the LTV 
framework still suitable for the analysis, we are looking forward 
for the answers. 



x{k+\) = Gk+ix{k) + Hk+xu{k) (18) 

Where G^+i =/-hr;t+i4+i> Hk+i = Tk+ib'-<'\ J^j^^^ and / 
are the transpose of the adjacency matrix of the {k+ l)th graph 
and the identity matrix, respectively. Define 



W, = [Gt--- G2Hu- ■ ■ ,GtHt-uHt] (19) 
and the final state is written as: 

x(r) =[xi(T),---,XN(T)] 

= [Gr---Gi]-[xi(0),---,xw(0)]' (20) 

+ Wciu(0)Mi),---MT-i)] 

If there exists a sequence of inputs denoted as 
[m(0),k(1), •••,M(r-l)]' such that [xi{T), ■ ■ ■ ,XNiT)]' = 
[0, • • • ,0] in Eq. (20), then the temporal network is structurally 
controllable at time point to, i.e. rank{Wc) = N. Otherwise, we 
may split x(T) into two parts, written as: 

x(T) = [xi(T), ■ ■ ■ ,Xk{T),Xk + \{T), ■ ■ ■ ,xm{T)] 

= [Gt--- Gi1-[xi(0), • • • ,Xkm,Xk+m, ■ ■ ■ ,xw(0)]' (21) 
+ W,-[u{Q),u{\),---,u{T-V)\ 

and if there exists a sequence of inputs denoted as 

KO),w(i),---,M(r-i)]' such that [xi(r),---xi:(r)]' = [o---o]' 

in Eq. (21), then the k subspace of the network is structurally 
controllable at time point ^o^ which is equivalent to the condition 
rank{Wc) = k. Therefore, we define controlling centrality as 

Suip) =rank{Wc) (22) 



Metliods 

4.1 Notation 

The symbols used in the main text are summarized in Table 2. 

4.2 Controlling Centrality 

With a samphng interval properly chosen, we write Eq. (3) as 
follow: 



=^^+^x{k)+b^''^u{k) (17) 



Generally, T^+i^Tk, where 7*:+i = ii+i — 4 is the sampling 
interval. From Eq. (17), we get the recursive relationship of two 
neighboring states as: 



i.e. the maximum dimension of controllable subspace, as a 
measure of node o's ability to structurally control the network: if 
Sm(o) = ^> then node o alone can structurally control the whole 
network. Any value of Sm(p) less than N provides the maximum 
dimension of the subspace o can structurally control. 

4.3 Datasets 

We mainly investigate three temporal networks with three 
empirical data sets in this paper. The first data was collected 
during the ACM Hypertext 2009 conference, where the 'Socio- 
Patterns' project deployed the Live Social Semantics applications. 
The conference attendees volunteered to wear radio badges which 
monitored their face-to-face interactions and we name this data as 
'HT09'. The second is a random data set containing the daily 
dynamic contacts collei:ted during the art-science exhibition 
'INFECTIOUS: STAY AWAY' which took place at die Science 
Gallery in Dublin, Ireland, and we name it as 'SG-Infectious'. 
These two data are both available from the website of 'Socio- 
Patterns' [31] (http://www.sociopatterns.org). The third data set 
was collected from Fudan University during the 2009-2010 fall 
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semester (3 whole months), which is named as 'FudanWIFF 
[35,36,43]. In this data set, each student/teacher/visiting scholar 
has a unique account to access the Campus WiFi system, which 
automatically records the device' MAC addresse, the MAC 
address of the accessed WiFi access point (APs), and the 
connecting (disconnecting) time as well. Table 3 summaries some 
characteristics of the aforementioned three empirical datasets. 
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